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Internal Energy of the Potts model on the Triangular Lattice with Two- 
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We calculate the internal energy of the Potts model on the triangular lattice with two- 
and three-body interactions at the transition point satisfying certain conditions for coupling 
constants. The method is a duality transformation. Therefore we have to make assumptions 
on uniqueness of the transition point and that the transition is of second order. These 
assumptions have been verified to hold by numerical simulations for q = 2, 3 and 4, and our 
results for the internal energy are expected to be exact in these cases. 

KEYWORDS: Potts model, exact energy, triangular lattice, duality transformation, transition 
point 

1. Introduction 

Exact solutions have exerted strong influences on the development of statistical physics of 
equilibrium systems, the spin models in particular. Among useful techniques for this purpose, 
one is the duality transformation, which applies mainly to two-dimensional cases. Ever since 
the first formulation of duality by Kramers and Wannier, 1 the method has been used to predict 
the transition points of the Ising model on the square lattice, 1 the Potts model on the square, 2,3 
triangular, and hexagonal lattices. 4 The duality transformation makes it possible to obtain 
the location of the transition point if it is unique. We can also calculate the internal energy 
at the transition point for a self-dual model from the derivative of the duality relation. 3,5 

In the presented paper, we apply the duality transformation to the Potts model with two- 
and three-body interactions on the triangular lattice and calculate the internal energy at the 
transition point. The Potts model with two- and three-body interactions on the triangular 
lattice was first introduced by Kim and Joseph. 4 They obtained this model as a result of 
an application of the duality and star-triangle transformations to the Potts model only with 
two-body interactions on the triangular lattice. Later, Baxter, Temperly and Ashley showed 
algebraically that the Potts model with two- and three-body interactions is self-dual, 6 and so 
also did Wu and Lin graphically. 7 It was also shown that the Potts model with only three- 
body interactions on the triangular lattice is self-dual, 6 though the Potts model with only 
two-body interactions is not. 

Baxter, Temperley and Ashley obtained the exact solution for the free energy and the 
internal energy of the Potts model at the transition point with two-body interactions on 
the triangular lattice. 6 They used the relationship between this Potts model and the soluble 
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Kelland model. However, it was difficult to analyze the Potts model with both two- and 
three-body interactions whereas the model is amenable to exact analysis by duality if it has 
only three-body interactions. 6 The Potts model with two- and three-body interactions has not 
been treated by such considerations, because of the complexity of the relationship between 
the original model and its dual although the duality for this model has been well known. 

Our goal is to derive the exact solution of the internal energy at the transition point for 
the Potts model with two- and three-body interactions using the duality transformation. We 
find that the duality transformation gives us the value of the internal energy only on some 
special points, which lie on the line of transition points. We also discuss the mechanism that 
we can calculate the internal energy by the duality transformation. 

The outline of the presented paper is as follows. Section 2 gives a brief introduction to the 
Potts model with two- and three-body interactions on the triangular lattice and its duality 
transformation. Here we will find the condition that we can calculate the internal energy by 
use of the duality transformation. In §3, we numerically verify the results obtained in §2. The 
last section is devoted to discussions. 

2. Energy of the Potts model with two- and three-body interactions 

We introduce the g-state Potts model with two- and three-body interactions on the trian- 
gular lattice and write its duality relations, following Ref. 8). The Hamiltonian is 

h(j 2 , Jz) = ^a(j 2 , j 3 ) = -J2^2Y, 6 (&> - J 3 E 6 Ofc> . (!) 

A A i^j A 

where {4>i} denotes the spin variables (fa = 0, 1, • • • , q — 1). The summation in the second 
expression is over all up-pointing triangles as shown in Fig. 1. Three-body interactions exist 
only on the up-pointing triangles. Therefore our model is slightly different from the model of 
Schick and Griffiths, 9 which has three-body interactions on all up-pointing and down-pointing 
triangles. The summation in the third expression with subscript i ^ j is over bonds (1,2), 
(2,3) and (3,1) around each up-pointing triangle as shown in Fig. 1. The strengths of the 
two- and three-body interactions are denoted as J 2 and J3, and 5 is the Kronecker's symbol. 
The partition function is written as 



Z(K 2 ,K 3 ) = £ exp {-(3H( J 2 , J 3 )} = £ ]J exp < 
{<Pi} {&} a 

(2) 

where K 2 = f3J 2 and K3 = f3J 3 . This model is self-dual and the duality relation is known to 

be 6 " 8 

Z(K 2 ,K 3 )=(^J Z(K*,K* 3 ), (3) 

where and are the dual coupling constants and N is the number of sites, which is the 
same as the number of up-pointing triangles. The relations between the coupling constants 
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Fig. 1. The Potts model with two- and three-body interactions on the triangular lattice. Each of the 
spin variables interacts with nearest-neighboring sites (two-body interactions). We also consider 
three-body interactions among spins on all up-pointing triangles as shown shaded. 



(K 2 ,K 3 ) and (K^K^) are 

v = — , y =—, (4 

y v 

where v = e K2 - 1 and y = e K[i+3K2 - 3e K ' 2 + 2 and similarly for v* and y* using K\ and K\. 
The fixed-point condition, K 2 = if| and K 3 = K%, is equivalent to y = q, which describes a 
line in the (K2,Ks) plane. It is expected that the system undergoes a phase transition across 
the line y = q. Indeed Wu and Zia argued that the line y = q represents the critical surface if 
J3 + 3J 2 > and J 3 + 2 J 2 > 0. 10 

Next let us calculate the internal energy on the line y = q, using the duality relation of the 
partition function (3). In the case of a simple single-variable duality Z{K) = c(K) N Z(K*), 
the logarithmic derivative of both sides readily leads to the energy at the fixed-point, 

1 dc 



e(K c ) = — 



s dK*\ dK 



(5) 

K=K C 



dK 

where e(K) is the energy per site, under the assumption of continuity e(K c + 0) = e(K c — 0). 
We therefore try to apply the same procedure to the present duality relation (3) with two 
variables. 

The logarithmic derivative of the partition function with respect to (3 gives the internal 
energy per site as follows; 

\ A #3 1 K 2 ,K 3 \ A I K 2 ,K 3 

= J 2 U 2 {K 2 ,K 3 ) + J 3 U 3 {K 2 ,K 3 ), (6) 

where the brackets (• • • )k 2 ,k 3 denote the thermal average. The values of U 2 (K 2 , K 3 ) and 
U 3 (K 2 , K 3 ) on the line y = q satisfy the following equations, which are obtained by taking 
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logarithmic derivatives of eq. (3) by K 2 and K3 and assuming continuity (i. e. second-order 
transition) of these derivatives across the line y = q, 



dig 

dKo 



Uo 



y=q. 



dK 3 



u* = - 



y=q 



'diq 



u 2 + l\ 



y=q 



diq 



y=i. 



dy 
dK 2 

dy 



y=q 



y=q 



(7) 



(8) 



It turns out that these two equations are not independent of each other as proved in Appendix. 
Because of this property, we cannot calculate both U 2 and U3 from eqs. (7) and (8), which 
makes it difficult in general to evaluate eq. (6). Nevertheless, if the left-hand side of eq. (7) 
is proportional to the final expression of eq. (6), eq. (7) gives the internal energy apart from 
the proportionality constant (which can be calculated from the values of coefficients U 2 and 
Us in eq. (7)). We therefore require that the ratio of coefficients of U 2 and U3 in eq. (7), 
{1 - (dK*/dK 2 )\ y=q }/{-(dK*/dK 2 )\ y=q }, equal to J 2 /J 3 : 



J A 1 



diq 
Wo 



-J 2 



y=i. 



dig 
dKo 



(9) 



y=q 



We can write this condition explicitly in terms of K 2 and K3, using the explicit forms of the 
derivatives in eqs. (A-4) and (A-6) in Appendix, as follows, 



(J3 + 3J 2 



v c (6v c + 3g) 

?(l+«c) 



J2 



3v c + q + 1 



3(1 + v c ) - (3v c + q) 



6v c + 3q 



(10) 



where v r 



D K 2 



1 with K 2 being the value on the line y = q. This equation is reduced to 



{K 3 + 3K 2 )(3e K2 +q- 2)(e K * - 1) = K 2 (3e K2 + 2q - 3)e K \ 



(11) 



where K3 is the value on the line y = q. We can calculate the internal energy at the limited 
points on the line y = q determined by eq. (11). 

Let us emphasize once again that the condition (11) for the internal energy to be calculable 
has been obtained under the assumption of uniqueness and continuity of the transition. The 
former uniqueness assumption was used in the identification of the fixed point y = q with the 
transition point as discussed in Ref. 10), and the latter continuity was necessary in deriving 
eqs. (7) and (8). 

3. Results and numerical verification 

Equation (11) with y = q has solutions for q = 2, 3 and 4 as well as for q > 69. The results 
for the former case are listed in Table I. The values of coupling constants in Table I satisfy the 
Wu-Zia criteria J3 + 3J 2 > 0, and J3 + 2J 2 > 0, which means that the existence of a transition 
in each case is very plausible. We have checked the assumption of uniqueness and continuity of 
the transition by numerical simulations for q = 2, 3 and 4. Systems with larger q are expected 
to have first-order transitions since the Potts model with two-body interactions on the square 
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q K2 K3 t energy 

2 0.375530 0.347552 0.925497 -3.19412 

3 0.291622 0.737724 2.52973 -3.95336 

4 -0.062177 1.75913 -28.2924 -15.5884 
4 -1.09861 4.39445 -4 -0.5 



Table I. Solutions of cq. (11) and y — q. Here r = K^jK-i. 



lattice has first order transition for q > 5. We thus have not carried out numerical simulations 
for the case of larger q. 

Most simulations have been performed for linear system size 100 and 10500 MCS per spin. 
The data have been averaged over 100 independent runs. The exception is the case with q = 4 
and r = K3/K2 = —4, in which we took 250 independent runs of 4500 MCS. 

The results are shown in Fig. 2 for q = 2, 3 and 4 (with r = -28.2924) and Fig. 3 
for q = 4 with r = —4. In Fig. 2 we observe that the internal energy shows no hysteresis, 
which suggests that the assumption of continuity of transition is valid. It is also seen that the 
duality predictions for the energy agree well with numerical data. In contrast, the case with 
q = 4 and r = — 4 seems to undergo a first-order transition as is seen in Fig. 3. Thus the 
solution given in Table I for q = 4 and r = — 4 would correspond to the average of e(K c + 0) 
and e(K c — 0), which is compatible with data in Fig. 3: In the case of a simple self duality 
Z(K) = c{K) N Z(K*), taking the limit K — > K c + after the logarithmic derivative leads to 
the average of e(K c + 0) and e{K c — 0) on the left-hand side of eq. (5). The same situation is 
expected for q = 4 and r = —4. 

4. Discussions 

We have calculated the internal energy of the Potts model on the triangular lattice with 
two- and three-body interactions satisfying certain conditions for coupling constants. These 
results have been obtained from duality under the assumptions on uniqueness of the transition 
point and that the transition is of second order. We checked these assumptions by Monte 
Carlo simulations. It has been shown that the assumptions are valid for q = 2, 3 and 4 
(r = —28.2924) but not for q = 4 (r = —4). In the former case it is expected that the duality 
predictions are exact. In the latter case our result gives the average energy above and below 
the transition point. 

It is possible to understand the condition (11) for the internal energy to be calculable from 
a little different point of view. In Fig. 4, r = K3/K2 and the dual r* = K^/K^ are plotted 
as functions of K2. As seen there, the two lines have the same slope (i.e. both are flat) at 
the fixed point. This feature dK^/dKZ, = dKj,jdK2 can be verified by direct manipulations of 
the conditions (7)-(ll). In this sense the solvability condition may be regarded as a condition 
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Fig. 2. The internal energy at the fixed points of duality as in Table. I except for q = 4, r = —4. The 
symbol 'inc' denotes the internal energy measured by increasing the temperature, 'dec' stands for 
the results by decreasing temperature, and 'dual' is for our duality results. 



for the dual model to be the same model as the original one, i.e., the same ratio of coupling 
constants between two- and three-body interactions. In other words the model is self-dual in 
a very strict sense only at the transition point for q = 2, 3 and 4 (r = —28.2924). 

It is interesting to see the case of q = 4, r = — 4 in Fig. 4: the model is self-dual even away 
from the transition point. This is a very interesting special feature of this case. 
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Fig. 3. The internal energy of the four-state Potts model with r = —4. The same symbols are used 
as in Fig. 2. 



Appendix: Evaluation of the determinant 

In this Appendix we show that eqs. (7) and (8) are not independent. The dual coupling 
constants K^ and K% are expressed by y and v as follows 



Kn 



K% 



log 1 + 



qv 

y 



-3K 2 * + log 1 + 



3qv + q z 

y 



where y = e K3+3I<2 — 3e K2 + 2 and v = e Ki — 1. On the line y = q 
v are seen to satisfy 

dy „ dy dv 



6v c + 3q, 



3v c + 1 + q, 



l + v c , 



the derivatives 

0v o, 



dK 2 dK 3 dK 2 dK 3 

where v c = e K2c — 1. From these relations, we can obtain the derivatives of if| and 
respect to K 2 and K% on the y = q line as follows 



dK 2 
dK* 



1 



dK 3 

diq 

~dK 2 



y=q 



y=q 



y=q 



y=q 



= -3 



_ v c (6v c + 3g) 
9(1 + v c ) 

v c (3v c + l + q) 
q(l + v c ) 

dK 2 



+ 



dig 

dK* 



y=q 



y=q 



3v c + q + 1 
(3v c + q) 



3(1 + v c ) - {3v c + q) 



6v c + 3q 



(Al) 

(A-2) 
of y and 

(A3) 



K% with 



(A4) 
(A5) 
(A-6) 
(A-7) 
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Fig. 4. The behaviors around the solutions of eq. (11). The horizontal solid line is the original model 
and the dashed line is for K^/K^. These two lines intersect at the fixed point (y = q) shown in 
black dot. 



Using these derivatives, the determinant of eqs. (7) and (8) is evaluated as 



1 - 



dKS 



dK 2 



1 - 



dKS 



y=Qj 



V=Qj 



dKo 



dKX 



y=q 



OK, 



v c (6v c + 3q) 



q(l + v c ) 

-v c (3v c + 1 + q) 
q(l+v c ) 



1 | 3 -^c(3t; c + l + g) _ (3v c + g) \ 
q(l + v c ) q J 



-3 1- 



v c (6v c + 3q) 
q(l + v c ) 



+ 



3v c + q + 1 



3(1 + v c ) - (3v c + q) 



6v c + 3q 

q 



0. 



Thus eqs. (7) and (8) are not independent. 
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